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We develop a general mapping from given kink or pulse shaped travelling-wave solutions including 
their velocity to the equations of motion on one-dimensional lattices which support these solutions. 
We apply this mapping - by definition an inverse method - to acoustic solitons in chains with 
nonlinear intersite interactions, to nonlinear Klein-Gordon chains, to reaction-diffusion equations 
and to discrete nonlinear Schrodinger systems. Potential functions can be found in at least a unique 
■ way provided the pulse shape is reflection symmetric and pulse and kink shapes are at least C 2 

functions. For kinks we discuss the relation of our results to the problem of a Peierls-Nabarro 
potential and continuous symmetries. We then generalize our method to higher dimensional lattices 
for reaction-diffusion systems. We find that increasing also the number of components easily allows 
for moving solutions. 
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• I. INTRODUCTION 

o 

I Finding exact travelling-wave (TW) solutions of nonlinear lattice systems has been a problem of growing interest 
in recent years. Apart from some integrable systems which support TW solutions (e.g., a little is known about 
non-integrable discrete systems. It appears to be difficult to prove the existence of such waves because one has to 
ON i deal with differential equations with advance and delay terms (on some general properties of these equations see 
H ) . The existence of acoustic (pulse) solitary waves as travelling- wave solutions in lattices with nonlinear intersite 
interactions has been proved in j|. However, no proof is available for other types of solitary waves, for instance, 
topological solitons in nonlinear Klein-Gordon (KG) lattices or other discrete kink-bearing systems (an exception is 
given in Q). Stationary breathers have been shown to be generic solutions for lattice systems (for a review and further 
references see [|) ) . Again the question of whether moving breathers on lattices exist is still not answered, although a 
number of approaches to the subject are known fij, 0, g, M], p| ). An exception is the case of the integrable 
Ablowitz-Ladik equation |fl2| . 

Here we approach the TW existence problem from the inverse side - we show that for a given TW profile, corre- 
sponding equations of motion can be generated, so that these equations of motion yield the chosen TW profile as a 
solution. This has been done first in ]l3j for the tanh-shaped kink and extended to reaction-diffusion-type systems in 
]l4j . However, in both cases the analysis was performed for a specific class of profiles, whereas we will approach this 
problem from a general point of view. This in turn will allow us to obtain general information about the properties 
of the TW solutions. 

The structure of the paper is as follows. In Section II we introduce the equations of motion. Section III is devoted 
to solutions of the nonlinear Klein-Gordon equation and of reaction-diffusion-type systems. In Section IV we study 
chains with nonlinear intersite interactions which admit acoustic (pulse) soliton solutions. We will refer to this type 
of lattices as to acoustic chains. In Section V we deal with discrete nonlinear Schrodinger-type (DNLS) equation. 
Section VI is devoted to the structural stability of solitary waves, and Section VII generalizes our method to higher 
space dimensions. Conclusions are given in Section VIII. 



II. THE EQUATIONS OF MOTION 



We consider a one-dimensional chain with lattice spacing equal to unity, which describes a system of interacting 
particles of unity mass. Such a system has a direct physical meaning and can describe, for example, simple quasi-one- 
dimensional molecular crystals. The interparticle interaction potential W n - n i (r) and the on-site potential V(u) are, 
in general, nonlinear functions: 
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(1) 



where u n is the displacement of the nth particle from its equilibrium position and m, n are integers. If the second 
derivative u in Eq. (Q) is replaced by the first derivative ii, we obtain a system of reaction-diffusion equations. 
Another system of interest is a generalized discrete nonlinear Schrodinger equation (DNLS) 

i<j>n + C(4> n+ x - 24> n + 4>n-l) 

+ F{\4> n \ 2 ){^ n -x + 0„+i) + G(|</>„| 2 )^ = (2) 

which appears in various fields. Here 4> n {t) is a complex- valued function and F and G are general nonlinear functions. 

We are not aware of any systematic approach which shows the existence or even obtains analytical expressions for 
TW solutions of the equations from above. Therefore we approach the problem from the opposite side. We formulate 
an inverse method of creating the potentials V or W or the pair of functions (F, G) for a given TW solution. 



III. SOLUTIONS OF THE NONLINEAR KLEIN-GORDON EQUATION 

For sake of simplicity let us consider the case of harmonic intersite interaction W n - n '(r) — [Cr 2 /2)8 n r n i±\y Then 
the equation of motion (Q) becomes the well-known nonlinear Klein-Gordon equation: 

u n = C(u n+1 - 2u n + Un-i) - V'(u n ). (3) 

First, let us study the dispersion law for small-amplitude waves oscillating around some minimum u m i n of the potential 
V(u). After linearising the on-site potential around the above-mentioned minimum the dispersion law can be written 
as follows 

u\ = 7 + 2C(l -cosg) (4) 

where q is a wave number and 7 = V"(u m i n ). The group velocity sq — duj q /dq attains its maximal value s max when 
dso/dq = 0: 



/ \/j 2 + 47C — 7 

Smax = y G - . (5) 

We are interested in TW solutions, i.e., solutions that propagate with a permanent shape and velocity: 

u n {t) = u(n — st) = u(z) , z — n — st , (6) 

where s is the velocity of the travelling wave. As a result, we obtain a differential equation with delay and advance 
terms: 

s 2 u"{z) = C[u{z + l)-2u{z)+u{z-l)]-V'[u{z)] . (7) 



A. Moving pulses 

First we consider solutions of a bell-shaped localized form, i.e., pulses. Given the profile u(z) and its velocity s, we 
can generate the on-site potential V. The function u{z) should satisfy the following conditions: 

• u(z — > ±00) — > 

• u(— z) = u(z) 

• u(z) is monotonic in [0, +oo[ 

• u(z) is analytic in [0, +oo[ 
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To show that the potential V can be generated in a unique way, we rewrite Eq. (^) in the following form: 

V'[u{z)] = f[u{z)] = -s 2 u"{z) 

+ C[u(z + 1) - 2<u(z) + u(z - 1)] . (8) 

Now we see that there exists a unique correspondence between the function u(z) and the force function f(u). Since 
u(z) is analytic, we can always rewrite u(z ± 1) in terms of u(z). Therefore for each u(z) with the conditions listed 
above the function f(u) can be uniquely defined. The potential V is then obtained by integrating f(u) once. This 
result does not change if we consider a more complicated interaction potential W which incorporates anharmonic 
terms and long range interactions. Thus, we showed here, that for a given interaction potential W, a given pulse 
profile u{z) which satisfies the above conditions and a given velocity of the pulse we always generate a unique on-site 
potential V which supports this TW profile as an exact solution of the equations of motion. 

What would happen if we loose the symmetry condition on u(z)7 Consider two values z\, z^ at which u(zi) = u{z2). 
Consequently, the argument of the force function is also the same. But the rhs of Eq. (||) will be different for z\ 
and zi in general, which implies that we obtain two different values for the force function at the same argument - 
a circumstance impossible for standard functions. Thus, we have to require the symmetry of u(z) which guarantees 
that the force function is defined in a unique way. This is also the reason why we can exclude the existence of more 
complicated pulse forms like anti-symmetric pulses, symmetric pulses with several maxima etc. 

Let us investigate some properties of V. Since for z — > oo u — > and u" — > 0, we find that /(0) = which was to 
be expected. To get more information about the dependence of f(u) for small u (which tells us about the stability of 
the TW solution) we need the leading order dependence of u on z for large z which is given by the ansatz of the TW 
profile. Let us assume that our ansatz yields an exponential decay of u(z) for large distances, i.e., 

u(z -> ±oo) - e-^ z l , fj. > . (9) 

After substituting Eq. ([)]) into Eq. (||) we obtain 

f(u) ~ -u\jj?s 2 - 2C(cosh^ - 1)] . (10) 

The slope of the force f(u) for small u changes its sign when C crosses the value C\ given by 

s 2 fi 2 

Cl = 2(cosh M -l)- 

This means that the potential V(u) has maxima at u — when C < C\ and minima if C > C\. Consequently, the 
asymptotic state u(z — > ±oo) = is a dynamically unstable one for C < C\ and a stable one for C > C\. There is 
another critical value [see again Eq. (|8|)] of C given by 

° 2 - 2[u(l)-«(0)] ■ (12) 

If C > C2, an additional extremum (maximum) in V(u) appears between u = and it(0). The possible scenarios 
Ci < C 2 and C\ > C 2 are shown in Fig. 1(a) and Fig. 1(b), respectively. The above statements about stability hold 
for any exponentially decaying pulse. If the decay is non-exponential, the on-site potential can become non-analytic 
at u = 0. For example, for a Gaussian tail 

u(z -> ±00) - e^ 12 ' (13) 

the on-site potential for small u gets dressed with logarithmic corrections, e.g. for fj, = 1 

f(u) = V'{u) ~ u [-s 2 (2 + 41nu) 
2C 



+ 2C cosh f 2yf- 

e V 



In u 



(14) 
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FIG. 1. The potential V(u) for a fixed pulse solution: (a) for Ci < C 2 and I: C < Ci < C 2 , 2: C\ < C < C 2j 3: 
Ci < C 2 < C; (b) for d > C 2 and 1: C < C 2 < C x , 2: C 2 < C < d, 3: C 2 < C x < C. The potential is obtained for positive 
u values and continued to negative values for sake of transparency. 



Note that the solution also exists in the "anti-continuum" limit C = 0. This seems to be surprising, since the 
oscillators are not interacting with each other. Still in this case we have a simple equation 



/[«(*)] 



/'(z) 



(15) 



It can be easily noticed from the bell-shaped form of u(z) that the function /(it) is anti-symmetric and has two 
zeroes, one of which is at u = 0. Thus, the potential V(u) has a maximum at u — and a minimum at some value 
< u < u(0) (see Fig. 2). The separatrix trajectory corresponds to the motion of each particle from the maximum 
of the potential V(u) to the right wall and back. This motion needs infinite time. 
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FIG. 2. The schematic representation of the pulse solution in the anti-continuum limit. 



It follows that it is possible to prepare the initial phases of all particles on this separatrix such that their uncorrelated 
motion resembles the motion of a pulse solution through the system. This solution is dynamically unstable because 
V"(0) < 0. 

Finally, we consider possible velocities for exponentially decaying pulses (^|). We want to check whether our pulses 
can be subsonic (s < s max ) or supersonic (s < s sma x)- Taking into account that [see Eq. (M)] 
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7 = C(cosh/i — 1) 



(16) 



we compare s with s r 
in this case 



Defining s = s ? 



s we find that the pulse can be both subsonic and supersonic since 



C 



cosh /z + 1 - v /(cosh^ + l) 2 -2(2 + /^ 2 ) 



(17) 



Consequently for fixed C and //, s has to be small enough to satisfy C > C\ [see (pi)], when s(Ci) > Sq. Thus for 
so < s < s(Ci) the solutions are supersonic, while for s < sq they are subsonic. 

Let us consider an explicit example of the sech-type profile. Suppose the profile is described by the function 



u{z) = uosech(fiz) 

where uq is the amplitude of the pulse and /i its inverse width. Using the expressions 

u"(z) — Mo/i 2 [sech(/^z) — 2sech 3 (/zz)] 
u 3 (z)~ 



u{z) - 2- 



cosh(/iz 

U{Z — 1) + U(Z + 1) = ZUn o ?f 



we reconstruct the on-site potential V(u): 

V(u) 



sinh (i + cosh (fJ,z) 

2 u{z) 

2u cosh jit — : — g 5 • 

«o + sinh u 2 (z) 



s 2 // 2 

r 2 2 i o/'Yi 2 i r* 4 



2« 2 



^ 9 cosh u 



sinh /x 



In 



sinh" /i 



(18) 



(19) 



(20) 



(21) 



It is easy to check that its shape will change with C exactly as described above. We can rewrite this potential in the 
following form: 



V{u) 



\k 2 u 2 + K2 n 9 2C u 4 + C\la? + u% In (a 2 + u 2 ) 



(22) 



where the parameters of the solution s and fi are given by 



(j, = arcsinh(uo/a) , s 



Vk 2 - 2C 
arcsinh(iio/a) 



(23) 



After proper rescaling constants C and uq can be eliminated and, consequently, we can reduce the number of system 
parameters to two. However, we are not able to rewrite the potential V(u) as a function independent from the 
parameters of the solution ([i"8"l ): s,Uo,fi. This means that we cannot answer so far whether the solution ( p"8| ) comes 
as a family of solutions of (123) or is a unique solution of the obtained equations of motion. 



B. Moving kinks 

Kinks or topological solitons are solutions which connect two minima of the on-site potential V(u). If V(u) has 
several equivalent minima, a countable infinite set of stationary (time-independent) kink solutions of the Klein-Gordon 
equation @ exists (in contrast to the space-continuous case, where the continuum groups of translation symmetry 
provides a smooth family of stationary kink solutions). Some of these solutions will be local minima of the total 
energy, and some will correspond to saddles. The question whether moving kinks as TW solutions (||) exist is still 
open. Some results (see, e.g., Q) suggest that kinks in discrete lattices experience a so-called Peierls-Nabarro 
barrier. One interpretation of this barrier is that it is the energy difference between stable stationary kinks and 



5 



unstable stationary kinks. Indeed it is clear, that to unpin a stable stationary kink, one needs at least this amount of 
energy. Another more sophisticated approach - the collective coordinate approach - is a projection technique which 
aims at accounting for the dynamics of a kink-like object boosted to move along the lattice. By coupling the kinks 
center of mass coordinates to phonons, one arrives at the result that a moving kink will radiate, loose kinetic energy, 
and finally be trapped (pinned) by the lattice. Here the barrier appears as a height of maxima in a potential which 
is used to describe the kink motion. 

Nevertheless, the analytical result of [|l3| suggests that it is possible to construct an exact moving kink solution. 
This result can be generalised for any profile u(z) that satisfies the following conditions: 

• u(z — > ±Oo) — > ±Uq 

• u{z) is monotonic in ] — 00, +00 [ 

• u{z) is analytic in ] — 00, +00 [ . 

Contrary to the case of pulse solutions, we do not need to require the function u(z) to have certain symmetries, so 
that we can restrict ourselves to monotonicity only. A non-symmetric u(z) profile will simply imply a non-symmetric 
function V(u). If the above-mentioned conditions are satisfied, we again can uniquely map the function u(z) onto the 
potential V(u). For simplicity, we renormalise the variable u(z) by the kink width u(z) — > u(z)/uq. 
If 

«(« ±00) ~ ± (l - e _ "W) , (24) 

we can perform the asymptotic analysis similar to the case of the pulse solution. If C < C\ [where C\ comes from Eq. 
(|Tl"|)] it follows V"(±l) < 0. Therefore the extremal points u = ±1 are maxima and, consequently, the asymptotic 
groundstates u — ±1 are unstable. The sign of V"(±l) changes at C = C\ so that for C > C\ these groundstates are 
stable. Another critical C2 value is given by 

_ a V"(0) 

° 2 _ 2[ U '(1)- U '(0)] ' 1 ' 

If C < Ci the state u = is a minimum and if C > C'2 it is a maximum and we have the standard double-well 
potential. For details see Fig. 3. 
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FIG. 3. The potential V(u) for a kink solution: 1: C < d < C 2 , 2: Ci < C < C 2 , and 3: Ci < C 2 < C; The potential is 
obtained for positive u values and continued to negative values for sake of transparency. 
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For 



Die, 



u(z) = tanh(/iz) (26) 
the explicit form of the potential V(u) was obtained in jl3| and |l4| and can be written as follows 

V(u) = (,s V - C)u 2 - ^fu 4 
C 

: — J— ln(cosh 2 fi-u 2 sinh 2 //) . (27) 

sinh /i 

The existence of such a solution does not imply the absence of Peierls-Nabarro barrier for stationary kink solutions. 
We demonstrate this for the particular case of the solution (p6|). We consider a chain with the on-site potential 
with parameters C = 1, s = 0.5, /i = 1. Then we calculate two stationary kink solutions, one of which is stat 
and the other one being a saddle in the energy landscape. Their energies are E\ = 0.991796 and E^ — 0.992153, 
respectively (measured relatively to the absolute energy minimum). The nonzero difference is the Peierls-Nabarro 
barrier. In fact, the existence of a Peierls-Nabarro barrier already follows from the stability of one of the stationary 
kink solutions, which implies that a finite amount of energy is needed to get out of the minimum. 

Furthermore, we now are in possession of a very effective method to generate equations of motion which support 
stationary kinks without a Peierls-Nabarro potential, i.e., where a stationary kink exists which can be placed anywhere 
in the chain. Such a system has a degenerated ground state under the constraint of existence of one kink. It can be 
easily generated by putting s — in the above inverse method, and that is also true for pulse solutions. Indeed, in 
the limit s — > we generate systems which support kinks (or pulses) which move with infinitesimally slow velocity. 
As the system is energy conserving and the kinetic energy becomes negligible in this limit, the groundstate becomes 
nearly degenerate, and a Goldstone mode appears for s — 0. We tested these predictions numerically and obtained 
excellent agreement. A particular example of a potential supporting kinks with zero Peierls-Nabarro potential is 

C 

V(u) = ~Cu 2 - - — — ln(cosh 2 n - u 2 sinh 2 y) , (28) 
sinh fi 

which was obtained for a kink solution u(z) = tanh(/Ltz) from ( p7|) by putting s — 0. It is interesting to note that Ward 
and Speight also proposed a scheme which generates systems supporting kink solutions with zero Peierls-Nabarro 
potential. This scheme uses Bogomolnyi's inequality ||l7| . The structure of the potential functions was fixed, but the 
difference operators were chosen in an appropriate way. This method then generates equations of motion which are 
rather hard to justify physically. This is not the case when using our inverse method scheme. 



C. Reaction-diffusion systems 

Let us also consider dissipative systems described by the discrete analogue of reaction-diffusion equation of the form 

il n = G{Un+l - 2u n + U„-l) + f(Un) . (29) 

The physical background of this equation differs from the above considered nonlinear Klein-Gordon chains but it also 
admits localized travelling-wave solutions. These systems are dissipative because we have a first order time derivative 
instead of the inertia term. Here the function f(u) can have different meanings, for example, an ion current for nerve 
fibers (see @). 

These systems lack time reversibility. The equation for TW solutions reads 

- su'(z) = C[u{z + 1) - 2u(z) + u{z - 1)] + f[u{z)] . (30) 

To generate a system for a moving pulse, let us consider a symmetric pulse with one maximum as described in the 
pulse section of the Klein-Gordon chains. Due to the first order derivative in Eq. (^0|) the lhs is anti-symmetric, 
while the rhs is symmetric. Consequently, we cannot define f(u) in a unique way. Any further complication of the 
symmetry of the pulse will not help either. Thus, we conclude that there exist no moving pulse solutions in Eq. (|3C|). 

However, reaction-diffusion systems ( pC| ) support kink solutions. The non-uniqueness problems disappear as long 
as the kink shape u(z) is a monotonous function. Due to the first order derivative in Eq. ( |3fj| ) a kink shape moving 
to the right with some given velocity will generate a function f r different from the function // generated by the same 
kink moving with the same velocity but to the left. 
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As an example, let us consider a profile u(z) = tanh(/iz). Performing the above-mentioned computations we obtain 

f(u) = -s^l-u 2 )-2C \ . + 2Cu (31) 

1 + (1 — u z ) sinn /i 

This result coincides with result of Bressloff and Rowlands (see Jl4| ) . 

Finally, let us show that we can also obtain moving pulses provided we increase the number of components per site. 
Indeed consider 

il n = C{u n+ i - 2u n + U n -i) + f(u n , V n ) , (32) 

v n = C{v n+ i - 2v n + v n -i) + g(u n ,v n ) . 

Assuming u n (t) — u(n — st) , v n (t) = v(n — st) we find 

-su'(z) = C[u{z + l)-2u(z)+u{z-l)] + f[u(z),v(z)} , (33) 
-sv'(z)=C[v(z + l)-2v(z)+v(z-l)]+g[u(z),v(z)} . (34) 

Let us choose a certain profile for u(z). Fixing a value of u — k we obtain a countable set of points such that 

u{zi) = k. Here i is an integer and counts all points. This defines a countable set of functions z\ u \u). Similarly 
we proceed with v(z). In order to solve the inverse problem, i.e., for given functions u(z) and v(z) we have only to 
require 

v(zl u) ) f v{zf) if j , u{zf ] ) ? u(z^) if i?j . (35) 

This is a weak condition satisfied by most choices of u and v. For instance we can even choose symmetric functions 
having just one maximum and decaying to zero at infinities. The only restriction would be to shift the centers of 
the two functions apart, e.g., v(z) = u(z — zq) for a given u(z). But also asymmetric functions with even several 
maxima are allowed. Also possible is a symmetric function for u(z) having one maximum and decaying at infinities, 
and an asymmetric function for v(z) having the same other properties - with the maxima of both functions coinciding 
(zq = 0). It is a tedious work to calculate examples, so in most cases it will be appropriate to obtain the functions 
/ and g numerically. The reason for the easy construction of two-component moving pulses is that we introduce two 
functions /, g of two variables, but determine them only on a line in their phase space {u, v}. That means that we do 
not completely define these functions. 

Adding a third component to the problem clearly further relaxes the conditions on the pulse forms. The existence 
of two component pulses is partially known for space continuous systems [ fL9[ . Note that our inverse method works 
as well in the space continuous case, i.e., where differences of the form u n +i — 2u n + u n — i are replaced by second 
derivatives. 

Why not doing the same for conservative systems? Then the functions /, g will be the components of the gradient 
of some generating function (e.g., a potential). Thus, we have to impose this gradient condition, which will restrict 
the choice of functions u.v. 



IV. ACOUSTIC CHAINS 



Now let us study systems which support acoustic (pulse) solitary waves. In these systems the on-site potential is 
absent [V(u) = 0] and the solitary waves appear due to the nonlinearity of the interaction potential W(r). First of 
all, we introduce the relative displacements, r n — u n +i — u n . In these terms the equations of motion take the form 

f n = W'{r n +i) - 2W'(r n ) + W'(r n -i). (36) 

For TW solutions r n (t) = r(n — st) = r(z) one can write 

W'[r{z + 1)] = sV'(z) + 2W'[r(z)\ - W'[r(z - 1)]. (37) 

As shown in ||, in such a lattice localized bell-shaped travelling- waves solutions can exist if W(r) has a hard anhar- 
monicity in the region r < 0. Note that the acoustic solitons correspond to a localized contraction of the chain and 
therefore the function r{z) should be completely negative. 

It is evident by following the above line of argumentation for Klein-Gordon chains, that the pulse r(z) has to be 
symmetric and must have only one maximum. Any deviation from this leads to a non-uniqueness in the definition 
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of the potential. This implies that acoustic chains admit at the best moving kink solutions (non-topological) in the 
original variables u n (t), and even these kink solutions have to have reflection symmetry. 

Suppose the function r(z) satisfies the conditions for the pulses given in Subsec. IIIA. Then due to the symmetry 
of r(z), 



W'(n) = W'(r ) + -sV'(O) 



(38) 



where tq = r(0) , r*i = ^(l). In order to find the unknown function W'(r) we have to solve the initial value problem 
( |37| ) where the initial condition should be a the function W' (r) on the interval ro < r < r\. If W'(r) is defined on this 
interval, using Eqs. (^) and (|38|), one can construct the function W'(r) for r > t\ which will be uniquely defined 
(see Fig. 4) for each W'(r). This means that we can choose W'(r) in [r ,ri] arbitrarily. 




r r l o 

r 

FIG. 4. The schematic representation of the creation of the potential W(r) for acoustic chains. Different initial value choices 
will generate different potentials for one and the same TW profile. 

Therefore we find that for given {r(z), s} a countable infinite dimension of the space of solutions W(r) £ C 2 exists. 
Each function from this set supports one and the same r(z) as an exact solution with one and the same velocity. 
However, the function constructed in this way from an arbitrary initial value will be non-analytic in general (it is easy 
to show that all functions will be twice differentiable). 

To avoid the problem of generating non-analytic potential functions, we found another way of constructing the 
potential. Suppose r(z) is a pulse. Then W'[r(z)] = p(z) is also a pulse. Let us rewrite Eq. ( p7| ) as 



r"(z) = 



p(z + 1) - 2p(z) + p(z - I) 



(39) 



Now instead of defining r(z) we define p(z) (symmetric bell-shaped pulse). If this function is analytic, the rhs of 
Eq. ( |39| ) is also analytic. Then simply integrating this rhs twice, we find an analytic function r(z) which is also of a 
bell-shaped form provided p(z) decays for large z faster than 1/z 2 . Having r as a function of p (on a half axis z) we 
can invert this dependence and consider p as a function of r. This gives us the force function W'(r). Notice that by 
that we can avoid generating non-analytic potentials. 

Finally, let us take a look at the asymptotic behavior of the solitary solution. Suppose p(z) ~ e~ M2 for z — > oo and 
W'(r) ~ Cr for r — ► 0. Linearising Eq. (6), we obtain 



s 

C 7 



cosh p — 1 

J 2 



> 1 



(40) 



Thus, each moving acoustic soliton is supersonic. 

In the following example we will illustrate how to construct the interaction potential from a given solution profile. 
Suppose p(z) = pq sech 4 (pz). We substitute it into the rhs of Eq. (p3) and integrate it twice. Using the formula 
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sech 4 (fiz)dzdz = — 



1 



6/i 2 cosh (fiz) 



2 

+ — — ln[cosh(^z)] + K\Z + K 2 



(41) 



we calculate the function r(z). In order to satisfy the boundary conditions, we put K\ = 0, K 2 = 0. As a result, we 
obtain 



Po 



3s 2 /i 2 



ln£(z) 



!p(z) 
Po 



2(sinhV- 1) 
TO 



, £(*) = ! + 



'^sinhV 
Po 



(42) 



We cannot express the force function W'(t') explicitly from the above formula, but we have the inverse relation 
r(z) = D[p(z)] where D is a function inverse to W' . The potential W(r) can be calculated numerically. 



V. BREATHERS OF THE DNLS-TYPE EQUATIONS 



Here we study a general nonlinear chain governed by Eq. (^|) . It is already known that these systems have standing 
breather solutions (see, e.g., p0|). The standing breather is defined as a spatially localized solution which is periodic 
in time. The general breather solution with frequency cj, velocity s and wavenumber q can be chosen to be 



<j> n {t) = $(n — st, fit) , $(x — > ±oo, y) — > , 
$(af,!/) = $(x,y + 27r) . 

Periodicity of <fr(x,y) in y allows to expand it into a Fourier series: 



(43) 



(44) 



For DNLS systems these solutions may have only one nonzero Fourier harmonic w.r.t. time. Since the DNLS equation 
has a gauge symmetry (f) n (t) — > e luJt 4> n (t), we can actually always transform a breather solution into a stationary pulse 
solution (note that this is not possible for breather solutions of e.g., Klein-Gordon or acoustic chains). As a result 
Eq. (0) can be rewritten as 



iijj n + K1p„ + C(lp n+1 + Ipn-l) + F(\ 1p n ?) (V'n-l + V'n+l) 

+ G(|^„| 2 )V„ = 



(45) 



where k = — 2C. Then we can define ijj n (t) — R(z) + il[z) , z — n — st , where the functions R and / are real. 
Separating real and imaginary parts of Eq. (|4q ) we obtain the unknown nonlinear functions F and G expressed in 
terms of breather envelope R and / and breather parameters 



F[R 2 (z) + I 2 (z)} = -C 
G[R 2 {z)+I 2 {z)] = — k - 



[R 2 (z) + I 2 (z)Y 



2 (R(z)[I(z - 1) + I(z + 1)] - I(z)[R(z - 1) + R[z + 1)]) ' 
[R\z) + P{z)]> 



R(z) 



I'(z)- 



2 (k*) - aigaffi ^) 



(46) 
(47) 



We are interested in symmetric profiles, i.e., w(z) = \/ ' R?(z) + I 2 (z) — u{—z) in order to ensure single value 
properties of the functions F, G. It is easy see that the profile is symmetric and bell-shaped if, e.g., both I(z) and 
R(z) are symmetric or one of these functions is symmetric and another one is anti-symmetric. The lhs of Eq. (B^) is 
symmetric. If R and / are symmetric, the functions R(z— 1) + R(z+ 1) and I{z— 1) + I(z+ 1) are also symmetric and 
the derivative of R 2 (z) +I 2 (z) is anti-symmetric. Consequently, the rhs of this equation is anti-symmetric. Obviously, 
this is possible only for a trivial solution ip n (t) = 0. 

Clearly moving solutions seem to have more complicated internal structures, as known from the solutions of the 
Ablowitz-Ladik system [O]. These solutions can be represented in the following form: 



4>n{t) = u n {t)e l ( qn -^ = u{z)e i{ - qn -^ 



z = n — st 



(48) 
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where u(z) is a real envelope amplitude, s is breather velocity, g its wave number and ui its frequency. Substituting 
this ansatz into Eq. (Eh, we obtain two equations for real and complex parts: 



(lu - 2C)u{z) + cosq(C + F[u 2 {z)])[u{z - 1) + u(z + 1)] + G[u 2 {z)]u{z) = , 
-su'(z) + sing(C + F[u 2 {z)])[u{z + 1) - u(z - 1)] = . 



(49) 
(50) 



After straightforward calculations the unknown nonlinear functions F and G can be expressed in terms of u and the 
system parameters: 



F[u 2 {z)] = -C-- 



u'{z) 



sing u{z + 1) — u(z — 1) 
G[u 2 (z)} = -{u - 2C) - cos q(C + F[u 2 (z)}) 



u(z — 1) + u(z + 1) 
u(z) 



-{u-2C) 



u'(z) u(z - 1) + u(z + 1) 



(51) 



(52) 



tang u(z) u(z + 1) — u(z — 1) 
If the envelope function u(z) satisfies the conditions 

• u{z — > ±oo) — > 

• u(— z) — u(z) 

• u(z) is monotonic in [0, +oo[ 

• u(z) is analytic in [0, +oo[, 

we can postulate again (similarly to Sections III and IV) that for any envelope u(z) defined as above and the set of 
parameters (s, q, u>) one can uniquely define the nonlinearity for the equation given by functions F and G. 



A. Examples 

Let us consider the particular case when u(z) = uosech 2 (fj,z) . Substituting this expression into Eqs. ( ]5l| ) and (p^), 
we obtain functions F and G: 



F(u 2 ) = -C 



sing sinh2/i 

2, . ' 2s l' 



sinh a 
1 H -u 



"o 



G(u J ) = -(w-2C7) 



tang sinh2/i 



, sinh u 

x cosh2/i it 

u 



(53) 



(54) 



J. T/ie Ablowitz- Ladik equation 



Let us look at the particular case when G = 0. In this case we have only one unknown nonlinear function, F. After 
simplifying the ansatz (|5l|)-(|5^) we obtain 



u'{z) (2C — w)tang u(z + 1) — u(z — 1) 



u(z) 

F[u 2 (z)} =-C + 



s u(z — 1) + u[z + 1) 

(2C-cj)u(z) 



cos g[w(z — 1) + u(z + 1)] 
In the particular case u(z) = UQsech(fxz) we obtain the quadratic function 



F{u 2 ) = -C + 



2C ~uj 



2cosg cosh^i 



sinh 



M 2 
— U 



(55) 
(56) 

(57) 
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We assume 



2C-uj 2C-u sinh 2 /i A 



= C , : 5-^ = - , A > . (58) 

2 cos q cosh /i 2 cos g cosh fi Uq 2 

Eq. (H) yields 

2C-w = SM - . (59) 

tan q tanh /i 

We can rewrite these equations in more common way, expressing the parameters of the solution s,u> and uo through 
q and [i: 

uo = y — — sinh/U , lu = 2C[1 — cosh/i cosq] , s = 2C— — sing . (60) 

V A /i 

This corresponds to the well-known integrable Ablowitz-Ladik equation. 



2. DNLS with local nonlinearity 
Now let us look at a well-known equation of the DNLS family. In the case of F(u 2 ) = 0, we have 

i4> n + C(4> n+1 - 24> n + ct> n -i) + G (|0„| 2 ) <t> n = . (61) 
We substitute the ansatz ( fl8| ) and consider Eq. (^) that for this particular case takes the following form: 

u'(z) = a[u(z + 1) - u(z - 1)1 , a = Cbmq , (62) 

s 

The absence of a to be defined function in this equation, contrary to the previous examples, makes this equation an 
equation for the pulse shape. Let us show that a pulse shaped function u(z) can not satisfy Eq. (|62|). Suppose first 
that our solution u(z) is periodic with some large period L. In this case we can expand the solution into Fourier series 

+°° / \ 

u(z) = u m exp(im—z\ . (63) 

m— — oo ^ ^ 

Substituting this expansion into Eq. ( |62] ) we obtain the algebraic equation 

2?r „ . / 2vr\ 
m— =2asmlm— 1 (64) 

where m is unknown integer. The equation x = 2a sin x always has always a finite number of roots for any non-zero 
a. Since m is integer, we can actually solve ( [m]) only for some specific values of a. This does not depend on L, so we 
can consider the limit L — > oo. A pulse solution u(z) would require an infinite number of nonzero harmonics in Eq. 
(|63|). Therefore it is impossible to satisfy Eq. ( |62] ) with a pulse shaped u{z). Consequently, Eq. ( |6l] ) does not admit 
moving breathers. 

It appears to be not possible to extend this method to systems like acoustic or KG chains. Note, that the spectrum of 
a DNLS-breather consists only of one frequency and therefore can be transformed into two differential delay/advance 
equations while acoustic/KG breathers will have an infinite number of harmonics and obviously cannot be rewritten 
as a countable number of retarded and advanced ODE's. 



VI. CONTINUATION OF MOVING SOLUTIONS FOR CONSERVATIVE SYSTEMS 



Let us discuss the question, whether a conservative system which allows for a certain moving solution and has been 
generated by our inverse method, has this solution as an isolated one, or as a part of a smooth family of solutions. 
In other words, we consider a given moving solution, generate the equations of motion, and search in the phase space 
vicinity of our solution for other moving solutions. That calls for a linear stability analysis of the phase space flow 
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around the given moving solution. Since our solution has some uniform asymptotic ground state (we assume that the 
parameters are such that the uniform asymptotic state is a ground state, see discussions above), we know that a part 
of the linear stability analysis spectrum will be given by just solving the eigenvalue problem of linearized fluctuations 
around the ground state. The eigenvectors will be plane waves e l ^ n ~ u i t ) (they will be deformed in the center of our 
moving solution) and their spectrum is given by some dispersion relation to q . Note that due to space discreteness ui q 
is periodic in q. Let us search for q values for which the plane wave can be cast into the form 



This is possible if 



qs = uj q 



(65) 



(66) 



Consider first the case of an acoustic chain. For small q we have ui q = vq with v < s (because all moving solutions 
will be supersonic, see discussions above). Consequently, there is only the trivial solution q = of Eq. (66), which 
simply corresponds to a shift of the center of mass of the acoustic chain. We can always work in the frame where the 
center of mass is resting at zero. Consequently, for acoustic chains we do not find plane waves which move with the 
same velocity as the original solution. That implies that there are no small perturbations of our solution which have 
plane wave asymptotics and yield again a moving solution. Thus, we conclude that if moving solutions are coming in 
families, all solutions on the family will be localized in space. 

Consider now the case of an optical chain (KG case). Since oj q= o ^ 0, we will always find at least one nonzero 
q- value (in general it will be a finite number of q values which tends to infinity as s — > 0) which solves Eq. ([66]). In that 
case we do know that there will exist perturbations with plane wave asymptotics which deform our original localized 
moving solution into a partly delocalized one, which is now in addition characterized by a nonzero amplitude in the 
asymptotics (i.e., the kink has oscillating tails at \z\ — > oo). Some numerical results that confirm the above-mentioned 
arguments are given in [ pl| . 

To further analyse this, we performed a numerical continuation of a moving solution using the pseudo-spectral 
method which is essentially a Newton method in q space (see, e.g., ^2|, We chose the moving kink solution 

( p6| ) which yields the potential ( p7| ) . The numerical method traces the phase space of the system for moving solutions 
nearby the starting one with slightly changed velocities and shapes. 



0.005 



0.025 




-0.025 



-0.005 



0.801 0.802 0.803 0.804 0.805 



FIG. 5. Dependence of the amplitude in the asymptotics of a moving kink on the velocity for a given equation (see text). 



In Fig. 5 we show the dependence of the amplitude A of the asymptotics versus velocity s. We indeed find that our 
chosen solution ( p6| ) (with C = 1, fx = In (1 + \/2), s = l/[\/21n (1 + \/2)j ~ 0.802278 in this particular case) can be 
continued, but it thus gets dressed with plane wave asymptotics. The change of sign in the amplitude implies a phase 
change from to n. A slightly changed potential will exhibit similar solutions but with slightly shifted curves in Fig. 5 
. Thus, it follows that the moving solution with uniform asymptotics is structurally stable, i.e., has a similar solution 
with uniform asymptotics for slightly changed equations of motion. This follows from the fact that the crossing of 
the A(s) curve in Fig. 5 with A = is a generic intersection. 
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VII. MOVING PULSES IN HIGHER LATTICE DIMENSIONS 



So far we have been discussing moving pulses and kinks in one-dimensional lattices. In this chapter we will show 
that the inverse method can be easily generalized to higher space dimensions for reaction-diffusion equations, provided 
we take into account two or more components. 

Let us consider a two-dimensional lattice. The function u n (t) will now depend on two lattice indices u mtn (t). The 
differences (it„-i + u n +i — 2u n ) will now turn into some general discrete Laplacians D(u m ^ n ). Assuming a moving 
solution in the form 

u m ,n(t) = u(m - s x t, n - s y t) = u(x, y) (67) 

we arrive at the equations 

- s x u tX - s y u t y + D[u(x,y)] = f(u) . (68) 

Fixing a value of u we obtain a line in the {x, y} space and since the rhs should not change, the lhs should be constant 
on this line - a very restrictive condition. If we instead consider two components moving in the same directions with 
same velocities, the equations become 

- s x u tX - s y u >y + D[u(x,y)] = f(u,v) , (69) 
-s x v, x - s y v t y +D[v(x,y)] = g(u, v) . (70) 

Again fixing a value for u we obtain some line in {x, y}. If we consider a functions u, v decaying to zero at infinities, 
this line will be a closed loop. Let us assume that v is not constant on the loops of constant u. That helps, but still 
if we fix some point on the loop with some given value of v there will be a countable number of other points Pi on 
the loop where v takes the same value. Then the lhs's of Eqs. ( |69j) and ([ro|) have to be equal in these points. We can 
satisfy this condition by demanding two symmetries but only in the case when we have only two points p^. First our 
pulse functions u, v should be invariant under reflections at a line parallel to the direction of motion. This ensures 
that the first order derivatives on the lhs's of Eqs. ( |69| ) and ( |70| ) will be the same in all pi. To ensure invariance of 
the Laplacians in pi we only have to demand that the chosen direction of motion (defined by s x ,s y ) is parallel to 
a reflection symmetry line of the lattice. For instance, for a square lattice these will be only the major lattice axes 
and the diagonals. The initially assumed condition that v varies along the loop of constant u can be, e.g., easily 
satisfied by considering pulses which are symmetric under point reflections and whose symmetry centers are shifted 
along the line of motion. Note that contrary to the one-dimensional case the inverse method yields the functions /, g 
in a two-dimensional part of their phase space {it, v}. 

What if we add a third component? The conditions weaken again, similarly to the case of two components in 
the one-dimensional lattice. For instance, one can design pulses where two components are invariant under point 
reflections with centers shifted along the line of motion, and the third component will be off-centered from the line 
connecting the two first centers. 

Let us consider a three-dimensional lattice and two components. Fixing u we now obtain a closed surface in {x, y, z}. 
Requiring v to generally vary on this surface, we find that v will stay constant at least on loops embedded on the 
surface. Since the lattice is invariant only under discrete symmetries, we can not satisfy invariance of the lhs's of 
equations similar to Eqs. ( |69| ) and ( |70| ) on this loop. Consequently, there exist no moving two-component pulses 
in a three-dimensional lattice (and straightforwardly in any higher-dimensional lattice). This is in contrast to the 
space-continuous case, where space is invariant under continuous symmetries. Then we can satisfy the invariance of 
the lhs's along the loop trivially if both pulses are invariant under rotations around a line pointing in the direction of 
their motion. The initial condition that v is constant only on loops (not on the whole surface) is easily obtained by 
considering pulses it, v with shifted centers, just as in the two-dimensional case. 

Adding a third component w to the three-dimensional lattice case reduces the problem of constant v and w on a 
loop to that in a countable set of points pi. Still we need a symmetry to ensure invariance of lhs's in the points pi. 
This is easily achieved in the case when we have only two points pi by demanding two symmetries. First we need all 
three components u,v,w to be invariant under reflection at a plane which contains the direction of motion. Secondly 
this plane has to be parallel to any mirror reflection symmetry plane of the lattice. Again the initial condition of 
having just countable sets of points pt with coinciding values {it, v, w} can be achieved by shifting the centers of the 
three pulse components apart while staying on one line - the direction of motion. For instance, for a cubic lattice with 
lattice points at x = l,y = m, z = n and l,m,n integer reflection symmetry planes are {x,y,0}, {0,i/, z}, {x,0,z}, 
{x, ±x, z}, {x,±z, z}, {x,y,±x\ among possible others. Any vector embedded in these planes is an allowed moving 
direction. 
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VIII. CONCLUSIONS 



In this paper we have studied several types of nonlinear lattice systems. Contrary to most of papers on nonlinear 
lattices where authors try to find a solution (either analytically or numerically) of the given system, we approach the 
problem from the opposite side - we look for the system (in fact, for the interaction and/or on-site potentials) which 
admits some specific solution. 

We have studied kinks and pulses in the Klein-Gordon system, acoustic solitons in chains with nonlinear intersite 
interactions and discrete breathers in the nonlinear Schrodinger-type systems. In all these cases the method enables 
us to generate a unique on-site or interaction potential for a given pulse or kink and its velocity if this solution 
satisfies certain conditions. As a particular result, we have shown that the acoustic solitons are always supersonic. 
We also conclude that nonzero Peierls-Nabarro barrier does not prevent discrete kinks from propagating with constant 
velocities. In the case of discrete moving breathers in DNLS-type systems we create nonlinear terms in the equation 
(Q) for given envelope profile and breather frequency, wave number and velocity. 

Our method is equally well suited for dissipative systems. Systems of coupled reaction-diffusion equations do not 
possess one important property which is time reversibility and therefore despite being closely related to the Klein- 
Gordon type equations, do not have pulse travelling- wave solutions in the one-dimensional case for one component. 
We generalize the search for pulses to higher lattice dimensions and find that moving pulses can be easily obtained 
provided we also increase the number of components. 

All presented results can be easily extended to systems with longer range interactions, and to space continuous 
systems (i.e., to partial differential equations). Note that the continuum limit of the considered difference equations 
is easily recovered by choosing solitary wave profiles which vary slowly along the lattice. 
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